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Abstract 

Conditional elimination of degrees of freedom is shown to lead to an exact 
expression for the rate of turbulent energy dissipation in terms of a renor- 
malized viscosity and a correction. The correction is neglected on the basis 
of a previous hypothesis Q that there is a range of parameters for which a 
quasi-stochastic estimate is a good approximation to the exact conditional 
average. This hypothesis was tested by a perturbative calculation to second 
order in the local Reynolds number, and the Kolmogorov prefactor (taken as 
a measure of the renormalized dissipation rate) was found to reach a fixed 
point which was insensitive to initial values of the kinematic viscosity and to 

values of the spatial rescaling factor h in the range 0.4 < h < 0.8. 
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In the numerical simulation of fluid turbulence, as in other areas of computational 
physics, there is a practical requirement to reduce the number of degrees of freedom ex- 
plicitly simulated. However, turbulence has the noteworthy requirement that any such re- 
duction must maintain the rate at which energy is dissipated. In this letter we put forward 
a method of renormalizing the energy dissipation rate. 

It seems to be widely understood that any attempt to reduce the number of degrees of 
freedom in the theoretical description of fluid turbulence requires some form of conditional 
average, in which the retained modes are kept unaveraged ||. Yet this requirement has not 
been recognised in most attempts to apply the dynamical Renormalization Group (RG) al- 
gorithm to turbulence. Normally such methods rely instead on a band-filtered unconditional 
average. A critical appraisal of some of the leading approaches in this area will be found in 
the paper by Eyink ||. 

Originally, our own work, although introducing some features of the conditional average, 
also relied on the use of the band- filtered unconditional average j|] . Later it was recognized 
that a conditional average in turbulence can only be carried out as an approximation, and 
the two-field decomposition was introduced to separate out random and deterministic effects 
Recently, we have redefined the conditional average in the form of a limit, eliminating 
the need to separate into two fields Q. 

The development reported here is that one of the corrections (to mode elimination) in 
the momentum equation, vanishes identically in the energy equation (and hence does not 
contribute to energy transfer), while a second correction contributes to the energy spectrum 
but vanishes identically in the equation for the dissipation rate. 

We consider homogeneous, isotropic, incompressible, stationary turbulence, with dissi- 
pation rate e given by 



where the approximate equality defines the maximum wavenumber fc . The value of k is of 
the same order as the Kolmogorov dissipation wavenumber k^ 1 = (e/z/g) 1 / 4 , where i/q is the 
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(1) 



kinematic viscosity. 

In such turbulence the pair-correlation takes the form (u a (k,t)up(k' ', s)) = 
D a p(k)Q(k;t, s)5(k + k'), where D a p{k) = 5 a p — k a kp/k 2 , and the energy spectrum is 
related to the spectral density by E(k,t) = A7rk 2 Q(k]t,t). To complete the specification of 
our problem, we assume that energy is being injected into some low range of wavenumbers 
by a source term W(k), which satisfies 

rdkW(k) = e } (2) 

J 

for some k <C fcjg . This ensures stationarity. 

Next we introduce a version of the RG which leads to a renormalized dissipation rate 
equation. The Navier-Stokes equation (NSE) may be written in dimensionless form as: 

{di+iy (k)k 2 }u a (k,i) = 

= R (k )M Q(3l (k) J d 3 j up(j, i)u^(k - j, i), (3) 

on < k < ko — 1, where k = k/ko, i = t/r(ko), u a (k,t) = u a (k,t)/V(ko), r(k ) is an, as 
yet, undefined timescale, V(k ) is the r.m.s. value of a velocity mode with \k\ = k , defined 
for any k by 

V 2 (k) = ^(u a (k,t)u a (-k,t)), (4) 

Ro(ko) = T~(k )V(k )kQ is the local Reynolds number (see Batchelor 0, pl07) and M a ^(k) = 
{2i)~ 1 [kfsD ai {k) + k^Dap^k)]. It should also be noted that the local Reynolds number is 
indeed non-dimensional, since u a (k, t) has dimensions L 4 T _1 . The dynamical RG algorithm 
can now be stated as follows: 

1. Rescale all wavevectors on k\ (=hk ), where < h < 1, for example k' = k/ki, such 
that < k' < h~ l , and then average out the effects of the high wavenumber modes to 
obtain a dynamical equation for the modes on the interval < k' < k[ (= 1). 

2. Use this low- wavenumber NSE to obtain the low- wavenumber energy balance equation. 
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3. Integrate the energy balance equation with respect to k! up to the value k' = k[ (= 1) 
to derive an equation for the dissipation rate. 

4. Repeat these steps until the dissipation rate reaches a fixed point, at a new maximum 
wavenumber k' N = h N k' Q . 

The first step is to write the velocity fields as 

u±(k,t) = V(k 1 )^(k',t'), (5) 

where ip~(k',t') is defined on < k' < 1 and i/j+(k', t') is defined on 1 < k! < ft, -1 . The 
combined (low- A; and high-fc) equation of motion then takes the form: 

{d tl + v' (k')k' 2 }^(k',t') = 

= J R 1 (fc 1 )M± 7 (k') J d»j / ^(j / ,t / )v 7 (fc / -i / ,tO, (6) 

where R^fa) = r(/ci) V(/ci)A:^. Next we average out the effect of the high-A; modes, while 
leaving the low- A; modes unaffected. In general this will require the conditional projection of 
some functional !F[il> a ] on the which we denote by a subscript 'c', viz. (-) c . This should 
not be confused with the usual ensemble average, as denoted by (■). An important property 
of the conditional average is that the constraint is lifted by a further unconditional average 

§• 

Taking the low- A; equation, as given by (6), we conditionally average both sides, and 
decompose the right-hand side according to (5), to obtain: 

{d t , + v' (k>)k'^-(k') = 

= R l (k l )M^(k') J dY{(^-(j> 7 -(fc'- j')) c + 

+ 2(^(f)^(k' - j')) c + (^(j')^(k' - j')) c }. (7) 

We shall see presently that only the last term on the RHS of (7) contributes to the energy 
dissipation rate. In |I| we gave a method for the approximate calculation of this conditional 
average. We shall refer to this approximation as the quasi- stochastic estimate (QSE) of 
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the conditional average and denote it by (-)qse- Then, rearranging (7) and adding and 
subtracting quantities to leave it unaffected, we may write the low-wavenumber equation as: 

(d t , + v>(k')k'^-(k')- 

-R^M-^k') J d 3 / (^{j')^{k f -j')) QSE 
= R^M-^k') J d 3 j'^(j')4j-(k' - j') + 

+ S~(k'\k[), (8) 



where 



S-(k'\k> 1 ) = R 1 (k 1 )M-^(k')x 

x / d 3 j'{(^(f)^(k' - j')) c - ^(j')^-(k' - f) + 



Si 

+ 2(^(j')^(k' -j')) c + 



s 2 

+ (^(j')^(k' - f)) c - (^(j')^(k' - f)) QSE }. 

" v ' 

S3 

(9) 

We now wish to obtain the QSE, (tpt^j'jt^ip^k' — j',t'))QSE, which appears on the LHS 
of (8). Following the procedure in Jl|, we form an equation of motion for this quantity from 
(6). This is solved perturbatively as a power series in R(k\) and the band-filtered moments 
of the . As a result, the low-fc equation, after eliminating the first band of modes, takes 
the form 

{d t , + v' (k')k' 2 }i;-(k') -J ds'A af3 (k',t' - s')^{k',s') 
= R 1 {k 1 )M- M {k f ) J d 3 j" (^(j')^-(k' - j'))c 

+ S-(k'\k[), (10) 



where 



A aP {k',t'-s') = D a(3 (k')Rl(h) x 

x[AW(k',t' - s') + A^\k\t' - s')R 1 (k l ) + 

+ A^(k',t'-s')R 2 1 (k 1 ) + ...}. (11) 



The coefficients A^°\ A^ l \ A^ 2 \ . . . depend on the moments of ip + of order 2, 3, 4, . . . respec- 
tively. It should also be noted that the even-order coefficients are real and the odd-order 
are imaginary, since the expansion is effectively in powers of M a p 1 (k), which is imaginary. 

We now form the energy balance equation for the explicit scales k < ki, by multiplying 
each side of (10) through by tp~{—k',t') and averaging unconditionally. We then multiply 
through by appropriate factors to restore the original unsealed variables, in the process 
introducing the energy spectrum E(k). We also add W(k), as specified in (2), with the 
result: 



where A(k) = tr A a p{k) and T(k,t) = J djT(k,j, \k — j\;t) is the usual transfer spectrum, 
with wavenumbers in the interval < k, j, \k — j\ < k±. 

Lastly, we may form an equation for the rate at which energy is transferred through the 
modes of the system. Integrating (12) with respect to k, we obtain 



Note that the integral over the transfer term vanishes identically due to the antisymmetry 



At this stage, all three renormalized conservation equations (for momentum, energy and 
dissipation rate) are exact. Now consider the effect of the term S~, divided into three parts, 
as shown in (9), and begin with the energy equation. The conditional average behaves as 
a stochastic variable under a further unconditional average ||. Thus for the first term we 
have (schematically) 



(d t + 2^(k)k 2 )E(k) +2 ds A(k, t - s)E(k, s) 



W(k) +T(k) + 8nk 2 V 2 (k 1 )(S-(k\k 1 )4>-(-k)), 



(12) 




(13) 



of f(k,j, \ k — j\;t) under interchange of k and j (see 0, p85). 



<SWC>~ / a 3 j'{({^(j')^(k'-j')U-(-k')) 
-(^(j'^-(k'-f)^-(-k'))} = o, 



(14) 
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and so the contribution from Si vanishes identically in the energy equation. 

Now we turn to the dissipation equation: evidently we need only consider S 2 and S3. 
The first of these gives 

J Sk> ~ 

~ / d 3 k'J d 3 j" (^(f)^(k' - j')^-(-k')} = 0, (15) 

by antisymmetry under interchange of k' and j'. It should be noted that this property holds 
only because both wavenumbers are on the same interval. This is not the case regarding the 
contribution from S3, which is 

|d 3 A;'(S^->~ 

~ / d 3 k'J dVi^W-jUl-*')) - 

- (^(j')^(k' - j')^-(-k')) QSE }. (16) 

However, note that the two terms will cancel under any circumstances in which the QSE is 
a good model for the exact conditional average. 

Now, in order to perform an RG-style iteration, we truncate the expansion for A a p, as 
given by (11), at lowest order and rename 

Rl(k 1 )A ( -°\k', 0) = 5v (k')k' 2 . (17) 

Note that A^> is imaginary and therefore cannot contribute to the dissipation rate. This 
means that we effectively neglect terms of order Rf(ki) and higher. To this level of approx- 
imation, (10) can be written as 

(d t/ +u' 1 (k')k' 2 )^-(k f ) = S-(k f \k[) + 

+ R^M-^k') J d 3 j%(j')^(k' - j'), (18) 

for < k' < 1. The renormalized viscosity is given by 

v[(k') = v> (k') + H(k'), (19) 
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and the equation for the increment is 

6v>(k>) = RKh) lim Ao(ifc') [ (fc/2/2 / ) 2 + C7 fc/j> 

I'^h- 1 [ j' z — K J jJL 

+ 0(Rt(h)), (20) 

where /x is the cosine of the angle between fc' and j' and Ao(fc') is the two-field form of the 
increment |J, given by 

A„rn- 1 fdV W'WfrW (2i) 

k' 2 J vb(j')j' 2 + v' (i')l' 2 ' { } 

Here f = 1 < < ^ and L(fe',j') = -2M^(k')M+ e (j')D+(k' - j'). 

Equation (20) for the increment to the viscosity involves lim \k! — j'\ — > (see |T]]), and 
if this limit is evaluated by taking the Q + (\k' — as an expansion in Taylor series about 
k' = 1, then we make contact with the two-field version of McComb and Watt ||, with 
spectral density given by 

«'^M^){ ftll/S -T ftl4/S (''-^)}- < 22 > 

It should be noted however that the factor in square brackets in (20) is new. This 
arises because we were able to improve on the Markovian approximation used in the earlier 
calculations of McComb et al. Details of this analysis will be given in a subsequent 

paper. The equations for any iteration labelled n can be found by induction, and numerical 
calculation shows that the renormalized viscosity approaches a fixed point for some n = N, 
where in practice N = 5 or 6, for most values of the spatial rescaling factor. At the fixed 
point, (13) may be rearranged to give 

e = 2 [ kN dkv N (k)k 2 E{k) - 
Jo 

-Sir dkk 2 V 2 (k N )(S 3 ^-(-k)), (23) 

where (Ssipai - &)) i s shown schematically in equation (16). 

Let us now consider how to assess this work. We begin by noting that the renormalized 
'viscous term' in (18) is not an observable, even though it may be calculated to any order 
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using (11). This is because S~ varies from realization to realization of the explicit scales. 
In contrast, the renormalized viscosity in (12) would be (if we took the step given in (17)) 
an observable, as all terms in the equation have been averaged. However, the contribution 
from S2 is likely to prove important for energy transfer and, as this has been omitted from 
our calculation, we shall leave discussion of this to a fuller account, and concentrate on (23) 
for the dissipation rate, as we know that S2 does not contribute to this equation. 

In order to calculate the renormalized dissipation rate, we assume a power-law form for 
the spectrum E(k) = ae r k s . Then, the requirement that the renormalized viscosity (19) 
and its increment (20) scale in the same way, along with the conservation requirement of 
(23), yield r = 2/3 and s = —5/3, along with an expression for the Kolmogorov prefactor a 
(see equation (92) in reference 0). 

The theoretical prediction of a can be taken as a measure of the predicted dissipation 
rate, and in Figure 1 we show the result of such calculations, with a iterating to a fixed 
point for several different starting conditions at one value of the spatial rescaling factor. 
The fixed point corresponds to the upper end of the inertial range and the value of a at 
the fixed point agrees well with the result obtained from numerical simulations. Figure 2 
shows the fixed-point value of a for a range of spatial rescaling factors. It is of interest 
to note that the chain-dotted curve to the left depicts an earlier version of the theory, in 
which the conditional average was approximated by a band-filtered average ||. The dashed 
line shows the result of working out the limiting form of the viscosity (the stochastic part) 
with scale separation ||, while invoking a Markovian approximation. The continuous line 
shows the effect of including the square brackets in (20) and gives rise to a prediction of 
a = 1.62 ± 0.02 for 0.2 < rj < 0.6, where the bandwidth r\ = 1 — h. Incidentally, it is 
perhaps worth remarking that the limit rj — > (which one would expect in the microscopic 
case) does not exist for macroscopic turbulence. This is a consequence of an exact symmetry 
of the NSE: local energy transfer vanishes when the wavevector triad takes the form of an 
equilateral triangle. 

Lastly, there is the question: how good is our perturbation calculation? The expansion, 
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which we truncate, is in powers of A = R^(k), with integrals over moments of the ip + (where 
iV^lrms < 1)- With the maximum value A = 0.16, this is a small parameter, but possibly 
not small enough for the truncation to qualify as a rational approximation ||. Accordingly, 
we may have to rely on the properties of the moment expansion. Certainly, the next step is 
to work out the magnitude of the fourth-order term. This is the subject of current work. 
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tute. We also thank Prof. M.E. Cates, G. Fullerton, A. Hunter and A. Quinn for reading the 
manuscript and making numerous helpful comments. CJ acknowledges the financial support 
of the Engineering and Physical Sciences Research Council. 
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FIG. 1. The Kolmogorov prefactor a reaching a fixed point for a variety of starting viscosities 
vq. (For the case where the spatial rescaling factor h = 0.60 or the bandwidth rj = 0.40.) 
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FIG. 2. Variation of the Kolmogorov prefactor a with bandwidth r? or spatial rescaling factor h. 
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